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1.  INTRODUCTION 

T-junctions  in  rectangular  waveguides  play  a  very  important  role  in 
the  design  of  microwave  circuits,  such  as  multiplexers,  power  dividers, 
directional  couplers,  filters,  and  phase  shifters  in  modern  communi¬ 
cation  systems  [1—4].  Since  the  modeling  of  T-junctions  is  a  classical 
problem,  a  variety  of  purely  numerical  techniques  or  analytical  and  nu¬ 
merical  methods  have  been  developed  during  the  past  decade.  Among 
them,  the  mode-matching  method  is  a  typical  one  of  analytical  and 
numerical  approaches. 

A  wide  range  of  waveguide  junctions  has  a  configuration  in  which 
several  uniform  waveguide  sections  are  connected  through  a  cavity  re¬ 
gion.  The  mode-matching  method  uses  the  expansions  of  the  fields  in 
the  waveguide  sections  in  terms  of  their  normal  modes.  The  expanded 
fields  are  matched  to  those  of  the  cavity  region  to  obtain  the  modal 
scattering  matrix  of  the  discontinuity.  This  requires  a  resonant  mode 
expansion  of  the  cavity  fields  [4,  5]  or  the  use  of  dyadic  Green’s  func- 
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tion  in  the  cavity  region  [6] .  In  order  to  avoid  such  a  sophisticated  field 
analysis  of  the  cavity,  various  techniques  [1,  2,  7,  8]  based  on  equiv¬ 
alent  circuit  concepts  have  been  incorporated  in  the  mode-matching 
procedure.  The  cavity  problem  with  multiapertures  is  reduced  [7]  to  a 
superposition  of  simpler  cavity  problems  by  subsequently  shorting  all 
apertures  but  one,  to  which  the  field  expansions  by  the  normal  modes 
can  be  applied.  This  principle  has  been  used  to  analyze  asymmetric 
series  E-plane  T-junctions  [1]  and  the  magic  T-junction  [2]. 

A  mode-matching  method  [3]  similar  to  [7]  has  been  devised  to  ob¬ 
tain  the  generalized  admittance  matrix  in  closed  form  for  E-/H- plane 
three-  and  four-port  junctions.  Although  the  method  provides  an  ac¬ 
curate  full-wave  analysis  for  rectangular  waveguide  junctions,  it  needs 
a  larger  number  of  modes  with  the  optimal  number  ratio  for  TE  and 
TM  modes  in  each  waveguide  port. 

In  this  paper,  we  present  the  mode-matching  method  combined  with 
the  Fourier  transform  technique  for  analyzing  rigorously  E-/H-plane 
multiport  junctions  in  rectangular  waveguides.  The  method  has  been 
recently  applied  [10]  to  the  analysis  of  E-plane  T-junction,  H-plane  T- 
junction,  and  right-angle  corner  bend.  We  discuss  here  its  extension 
to  the  problem  of  E-/H-plane  cross  junctions.  In  this  approach,  the 
multiport  junctions  are  treated  as  a  problem  of  rectangular  apertures 
located  on  the  E-/H-planes  of  a  main  waveguide.  The  electric  and 
magnetic  fields  in  the  main  waveguide  expressed  in  terms  of  Fourier 
integrals  are  matched  on  the  apertures  to  those  of  arm  waveguides 
expanded  in  terms  of  normal  modes.  The  mode-matching  procedure 
is  performed  in  the  Fourier  transformed  domain.  This  yields  a  system 
of  linear  equations  in  closed  form  which  relates  the  modal  expansion 
coefficients  in  the  arm  waveguides  to  the  field  of  initial  excitation.  The 
main  advantage  of  the  method  is  that  the  scattering  parameters  are 
calculated  at  one  time  by  solving  matrix  equations  of  relatively  small 
dimensions  without  the  relative  convergence  problem. 

The  proposed  method  is  applied  to  the  analyses  of  the  Magic-T  junc¬ 
tion  and  an  asymmetric  E-/H-plane  three-port  junction.  It  is  shown 
that  the  convergence  of  the  numerical  solutions  is  very  fast.  The  re¬ 
sults  are  compared  with  available  numerical  and  experimental  data 
[2,3].  The  good  agreement  between  them  confirms  the  validity  of  the 
present  method. 
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2.  THEORY 

The  configuration  of  an  E-/H- plane  four-port  junction  is  shown 
schematically  in  Figure  1,  which  consists  of  three  rectangular  waveg¬ 
uides.  For  convenience,  the  waveguide  I  extending  from  the  port  4  to 
the  port  1  along  the  z  axis  is  referred  to  as  the  main  waveguide,  and 
the  waveguides  H  and  HI  comprising  the  ports  2  and  3,  respectively, 
are  referred  to  the  arm  waveguides.  The  cross  sectional  dimensions 
of  the  waveguides  I,  I  and  HI  are  2a  x  26 ,  2  ax  2 1  and  2  w  x  26 , 
respectively.  We  assume  that  the  TE\q  mode  expressed  by  the  longi¬ 
tudinal  section  TMX  waves  is  incident  from  the  right  side  end  of  the 
main  waveguide  I.  When  the  left  side  end  of  the  main  waveguide  I  is 
short-circuited  at  2  =  0,  the  four-port  junction  is  reduced  to  an  E-/H- 
plane  three-port  junction.  To  treat  the  problems  of  the  four-port  and 
three-port  junctions  by  the  same  formulation,  we  express  the  electric 
Hertzian  vector  of  the  incident  field  as  follows: 

n  fn  =  ^  sin  ai(y  +  a)  \ej0oZ  +  Te~jf)oZ]  (1) 

k0  L  J 

where  a\  =  vr / (2a)  ,  (3q  =  k q  —  a\  ,  ko  =  cu^/eo/io  is  the  wavenumber 
of  free  space,  T  =  0  for  the  original  four-port  junction,  and  T  =  —  1  for 
the  modified  three-port  junction.  The  first  term  in  the  right-hand  side 
of  (1)  denotes  the  original  incident  wave  with  a  unit  amplitude,  and  the 
second  term  with  T  =  —  1  gives  the  reflected  wave  caused  by  the  short- 
circuited  conductor  at  2  =  0.  The  scattered  fields  into  arm  waveguides 
H  and  HI  are  expanded  using  the  normal  modes  of  ( TEy  ,TMy )  and 
( TEX ,TMX )  of  each  waveguide,  respectively,  as  follows: 

OO  OO 

n!  =  §2  X]  XT  Amne~jkmn(-X~b )  cos  am(y  +  a)  sin  t„(z  -  d  +  t ) 

o  771=0  77=1 

(2) 

^  OO  OO 

nn  =  jy^7~  Bmne~jkmn(-X~b)  sin  am{y  +  a)  cos  tn(z  -d  +  t) 

(3) 

zv  OO  OO 

n&  =  §2  X]  ^2  sinmm(2  -  c  +  w)  cos  bn(x  +  6) 

o  777=1  77=0 


(4) 
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Cross  sectional  view  in  x=b  plane  Cross  sectional  view  in  y=a  plane 


Figure  1.  (a)  Schematic  of  an  E-/H- plane  four-port  junction  and  its 
cross-sectional  views  (b)  in  the  plane  x  =  b  and  (c)  in  the  plane  y  =  a. 

^  OO  OO 

nm  =  tX~  X  X  cos  wm{z  - c  +  w )  cos  bn(x  +  b ) 

o  o  m=Q  n= i 

(5) 

where  am  =  m7r/(2a)  ,  wm  =  mir/(2w)  ,  bn  =  nir/(2b)  ,  tn  =  nir/(2t)  , 

kmn  =  \lkl-a2m-  t2n  ,  gmn  =  \Jkl-w2m-  b„,  Z0  =  \J no /eQ  is  the 
intrinsic  impedance  of  free  space,  and  Amn,  Bmn,  Cmn  and  D.mn  are 
unknown  expansion  coefficients  of  the  transmitted  waves. 

The  Fourier  transform  technique  is  used  to  represent  the  scattered 
fields  in  the  main  waveguide  I.  The  technique  can  be  directly  applied 
to  the  problem  of  the  four-port  junction  with  the  main  waveguide  of 


Analysis  of  E-/H-plane  junctions 


277 


an  infinite  extent  in  the  z  direction.  For  the  problem  of  the  three- 
port  junction  with  the  short-circuited  conductor  at  z  =  0 ,  we  intro¬ 
duce  an  image  waveguide  structure  [10]  in  z  <  0  and  transform  the 
semi-infinite  waveguide  I  into  an  infinite  uniform  waveguide.  Using  a 
similar  analytical  procedure  as  described  in  [10],  the  electric  and  mag¬ 
netic  Hertzian  vectors  representing  the  scattered  waves  in  the  main 
waveguide  I  are  expressed  by  the  Fourier  integrals  as 

i  r°° 

Znkl  7^00 

yBl{Q  sm£„(x  +  b )  cos  au(y  +  a)^e~j^zd(  (6) 

1  f°°  f  - 

2^Z0L J^O^Kix  +  b) 

yB„(C)  cos£j,(x  +  b)  sin au{y  +  a)}e~j<’zd(  (7) 

where  ry,  =  \J  k^  —  —  £2,  £„  =  \Jk\  —  a2  —  £2  ,  and  an  infinitesi¬ 

mal  small  loss  has  been  assumed  in  the  wavenumber  ko  of  free  space 
which  is  finally  reduced  to  zero.  In  (6)  and  (7),  Af,(£)  and  A^(() 
represent  unknown  spectral  functions  for  the  scattered  fields  produced 
by  the  discontinuity  on  the  E-plane,  whereas  and  B%( Q  repre¬ 

sent  those  by  the  discontinuity  on  the  H-plane.  Note  that  the  electric 
fields  derived  from  (6)  and  (7)  satisfy  the  boundary  conditions  on  the 
conducting  walls  at  x  =  —  b  and  y  =  —a. 

The  tangential  components  of  electric  and  magnetic  fields  derived 
from  (l)-(7)  should  be  continuous  across  the  boundary  planes  x  =  b 
and  y  =  a  .  Taking  into  account  the  relations  of  symmetry  [10]  of  the 
fields  in  the  original  and  image  waveguides,  these  boundary  conditions 
may  be  expressed  as  follows: 


f  En,y(b,y,z) 

|y|  <  a,  \z  —  d\  <t 

Ei,y(b,y,z)  =  1  TEkty(b,y,-z) 

\y\  <  a,  \z  +  d\  <t 

(8) 

l  0 

otherwise 

f  E^z(b,y,z) 

\y\  <  a,\z  —  d\  <  t 

B\,z(b,  y,  z)  =  <  -TEk,g(b,y,-z) 

\y\  <  a,  \z  +  d,\  <  t 

(9) 

l  0 

otherwise 

f  EWtX{x,a,z) 

|x  <  6,  \z  —  c  <  w 

E\,x{x,  a,  z)  =  <  TEW:X(x,a,-z ) 

|x|  <  6,  \z  +  c  <  w 

(10) 

lo 

otherwise 

| xA%(Q  cos  bv(x  +  b)  sin  ry,(y  +  a) 


nr  =  £ 

i/=0 

-1 

oo 

n  f  =  E 

u=0 

4 
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{Ew,z(x,  a,  z )  |x|  <  b,  \z  —  c|  <  w 

— r'-Em.zC^j  a,  —  z)  |a;j  <  b,  \z  +  cj  <  w  (H) 

0  otherwise 


HTL,y(b,y,z) 

=  Hlty(b,  y,  z)  +  Hin,y(b,  y,  z) 

A 

a 

TT 

1 

a. 

A 

to 

Hj,z(b,y,z) 

=  HitZ(b,  y,  z)  +  HintZ(b,  y,  z) 

|y|  <  a,\z  —  d\  <t 

(13) 

Hw,x{x,a,z ) 

=  Hi.x(x,  a,  z)  +  H„hX(x,  a,  z) 

\x\  <  b,  \z  —  d \  <  w 

(14) 

Hw,z(x,a,z) 

=  HijZ(x,  a,  z)  +  Hin>z(x,  a,  z) 

|x|  <  b,\z  —  d\  <  w 

(15) 

where  Hin^x  to  Hmz  represent  the  components  of  magnetic  field  of 
the  incident  wave  deduced  from  (1).  The  boundary  conditions  for  the 
electric  fields  are  first  applied.  From  (6)  and  (2),  the  ^components  of 
the  electric  fields  on  the  plane  x  =  b  are  derived  as  follows: 


Ei,y(b,y,z) 

oo  /*oo  j  2  _ 2 

=  Y  V72  /  Bt(0  -  Y"  Sin(26^)  cos  a„(y  +  a)e~Ezd( 
u=0  z-KKq  J  _00  K q 

(16) 


(17) 

Equations  (16)  and  (17)  are  substituted  into  the  boundary  condition 
(8).  Then  (8)  is  integrated  from  y  =  —a  to  y  =  a  after  multiplying 
both  sides  by  the  trigonometric  functions  cos  aM(y  +  a)  ,  where  y  is 
nonnegative  integer.  The  Fourier  transform  of  the  resulting  expression 
is  calculated  with  respect  to  z-  coordinate.  This  leads  to  an  equation 
which  relates  the  unknown  spectral  function  E>®(£)  to  the  expansion 
coefficient  A/in  as 


Ei,y  (b,y,z)  = 


E  E  -4" 

m= 0 n=l 


avq  am 
k0 


cos  am(y  +  a)  sin  tn(z  —  d  +  t) 


bu  o  =  YA^ 

n= 1 


tn©n (d,  C) 

sin(2^6) 


(18) 
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where 


&n(d,t,  C) 


Un(d,t,C ) 


Un(d,t,C)  +  TUn(d,t,-C) 


eKd  [(-1  )■*&&  -  e-Kt' 

e-ti 


(19) 

(20) 


The  z-component  of  electric  fields  derived  from  (2),  (3), (6)  and  (7)  are 
substituted  into  the  boundary  condition  (9),  and  the  similar  procedure 
described  above  is  applied  to  obtain  the  following  relation: 


71=0 


k/j,nC®n(di  t ,  0 

Osin(2  ^b) 


(21) 


The  boundary  conditions  (10)  and  (11)  are  applied  in  a  similar  way. 
The  x-  and  z-conrponents  of  the  electric  fields  derived  from  (4)— (7) 
are  substituted  into  (10)  and  (11),  respectively.  For  this  case,  the 
trigonometric  functions  cos  b^x  +  b)  are  multiplied,  and  (10)  and  (11) 
are  integrated  from  x  =  —  b  to  x  =  b .  The  Fourier  transforms  of  the 
resulting  expressions  lead  to  the  equations  which  relate  the  spectral 
functions  Ae^(Q  and  A^{C)  to  the  expansion  coefficients  Cmfl  and 
Dmi, ,  respectively,  as  follows: 


OO 

^(0  =  E  cw 

m—  1 
oo 

A%)  =  E  D"w 

m=0 


Wm@  m  (c,w,Q 

sin(2  7/Ma) 

Qm.fiC®  m  (c,w,0 
Vn  sin(2 rj ^a) 


(22) 

(23) 


Next,  the  boundary  conditions  (12)— (15)  for  the  magnetic  fields  are 
applied.  From  (1)— (7),  the  ^components  of  the  magnetic  fields  on  the 
plane  x  =  b  are  given  by 


Hhy(b,y,z) 


'vVv  Ah 


k0 


4X0 


sin  rjv(y  +  a)e  ^zd( 


k20-C2 

ko 


4X0  cos(2&£„)  sin  a„(y  +  a)e  j<^zd( 


(24) 
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JPo  z  _  Ye~j0oz 


Hin,y(b,y,z)  =  —  -^-7—  sin  a\ (y  +  a) 

Zofco 

00  00  k2  -  a2 

Hn,y{b,y,z)  =  ^2  Bmn  0  ,  2m  sin  am(y  +  a)  cos  tn(z 


m= 1  n= 1 


(25) 

d  + 1) . 

(26) 


Equations  (24)-(26)  are  substituted  into  the  boundary  condition  (12), 
and  the  unknown  spectral  functions  A®(£)  ,  -A^(C)  1  and  B^(C)  are 
replaced  by  the  series  of  the  modal  expansion  coefficients  given  by 
(18)-(23).  Then  (12)  is  integrated  over  \y\  <  a  and  \z  —  d\  <  t 
after  multiplying  both  sides  by  the  trigonometric  functions  sin  ap(y  + 
a)  cos  tg(z  —  d  +  t) ,  where  p  and  q  are  nonnegative  integers.  This 
leads  to  a  set  of  linear  equations  for  the  expansion  coefficients  Bmn  , 
Cmn  ,  and  Dmn  as  follows: 


ko  ~  al 


T 2 Bpqt(l  +  5go) 

/Cn 

j(ko  R  T  y^  j(-l)p+napwm 

/  J  K'pn-E'pn-L  1  /  J  /  J  ^mri-L 2 


k2 

n= 0 


m=  1  n=0 


+  f;  f;  + ^Mneq(d,t,M 


m=  1  n=0 


7  I  7 

kq  a  ko 


(27) 


for  p  =  1,  2,  •  ■  • ,  <7  =  0,1,2,- 


with 


A 


/2 


(2@n{d,t,()Uq{d,t,-()  ^ 

27r/coCp  tan(2£p6) 
C2Qm(c,  w,  Q^g(d,  T  -C) 
2Trko(ri2  —  a2) 


(28) 

(29) 


where  5mn  denotes  the  Kronecker’s  delta.  The  integrals  in  (28)  and 
(29)  are  evaluated  in  closed  form  by  a  simple  residue-calculus  (See 
Appendix).  Applying  the  same  procedure  as  described  above  to  (13), 
we  have  another  set  of  linear  equations: 
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kpq  .  aptq 

ko  pq  kl  pq_ 


f(l  +  5po) 


(1  +  <5po)  'y  ^  jtntqApnls  j(l  +  <5po)  ^  ^ 


CLptqkpn 

~  1 


n=  1 


77=0 
oo  oo 


fco 


v  e  -  e  f>(- i)"+”5=^iDm„/2 


777 =1  77=0 

k0 


@q(d,t,Po) 
for  p  =  0, 1, 2,  • 


777=  1  n=0 


<7  =  1,2,  ••• 


/SQd 


(30) 


where 


h  = 


I 4  = 


ZPen(d,t,QUq(d,t,-  c) 


27r/co  tan(2£p&) 


dc 


poo  __2 


Vi@m(c,W,C)Uq(d,t,-t) 
2nk0{pl  -  a^) 


dC- 


(31) 

(32) 


Finally  we  apply  the  boundary  conditions  (14)  and  (15).  The  x-  and  z- 
components  of  the  magnetic  fields  on  the  plane  y  =  a  derived  from  (1)- 
(7)  are  substituted  into  (14)  and  (15).  In  this  case,  (14)  and  (15)  are 
multiplied  by  the  trigonometric  functions  sin  wp(z  —  c+ w)  cos  bq(x  +  b ) 
and  integrated  over  \x\  <  b  and  \z  —  c|  <  w  .  Then  we  have  other  two 
sets  of  equations  for  the  expansion  coefficients  as  follows: 


ko  ~  b2q 


wb(\  +  5po)Dpq  —  EE  j(-l)m+%inAmnJ2 

0  m=0n=l 

fcm.n&mfoi 
k0 


t2x 


+  -  jb  e  g"m((°  i^-Dmnjl 


m=  1  77=0 


777=0 


ko 


(33) 


for  p  =  0, 1, 2, 


9=1,2,' 
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wpbq  p  ,  -Jpq  ri 
k2  pq  +  fc0  ' 


pq 


wb(  1  +  Sqo)  —  —  EE  j(-l)m+HnwpAmnJ4 


m= 0  n= 0 


oo  oo  ,  _  _  oo 

+  EE  j(  1)’"+’  BmnJ2  +  1,(1  +  J0o)  E  jWm^p  Crnq  J 3 

m=l  n=0  m=  1 


+  6(1  +  <5go)  ^2  j 


■  Sm.nbqWp  2jbcL\Wp5qQ 


Dmq  J 1 


m= 0 


k0  k0 

for  p=  1,2,---,  g  =  0,1,2,- 


0p(c,w,/3o) 


(34) 


with 


Ji 


J2 


J4 


(2@m(c,  W,  C)Up(c,  W,  ~ C)^ 
2-Kkorjq  tan(2r]qa) 

(2@n{d,  t, QUpjc,  w,  -Q  ^ 
2^o(£m  -  &J) 

r)qQm{c,  w,  QUpjc ,  w,  -Q  ^ 
27rfco  tan(2?yga) 

£m®n(d,  t,  QUpjc,  w,  -C)  ^ 
2vrfc0(^  -  6g) 


(35) 

(36) 

(37) 

(38) 


where  Ji  to  J4  denote  the  integrals  characterizing  the  junction  in 
H-plane  and  are  the  counterpart  of  I\  to  I4  defined  by  (28),  (29), 
(31)  and  (32)  for  the  junction  in  E-plane.  The  integrals  in  (35)— (38) 
are  easily  evaluated  in  closed  form  by  a  residue-calculus,  though  the 
details  have  been  omitted. 

Equations  (27),  (30),  (33)  and  (34)  are  solved  to  obtain  the  un¬ 
known  expansion  coefficients  Amn  to  Dmn  for  the  scattered  waves 
into  waveguide  H  and  HI ,  after  truncating  the  modal  expansion  up 
to  m  =  n  =  M .  The  results  are  used  in  (18)  and  (21)-(23)  to  deter¬ 
mine  the  unknown  spectral  functions  -B)j(C)  to  A^(Q  .  Using  (18)  and 
(21) — (23),  the  Fourier  integrals  in  (6)  and  (7)  are  evaluated  in  closed 
form  by  a  simple  residue-calculus.  The  parameter  T  discriminating 
between  the  four-port  junction  and  the  three-port  junction  is  included 
in  the  function  Qn(d,t,()  defined  by  (19).  We  may  take  T  =  0  for 
the  four-port  junction  and  T  =  —  1  for  the  three-port  junction  with 
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a  short-circuited  conductor  at  z  =  0.  When  the  junction  is  excited 
by  TE\q  mode  incident  from  waveguide  I,  the  scattering  parameters 
«S'n,  52i,  531 ,  and  541  are  calculated  in  terms  of  expansion  coefficients 
rooted  by  (27),  (30),  (33)  and  (34)  as  follows: 


5n 


521 

531 


54i 


M 

E 

m= 1 


- jaiWm 0  m  (c,  w,  (3o) 


2a/30 


c  |  ^2  jklnPo@n{d,t,(30)  B 


n=0 


4bko 


In 


10 


I  a(fcg  —  w\ )  2 


C10 


(39) 

(40) 

(41) 


EhI  -jaiwmUm(c,  w,  -Po)  jkinp0Un(d,  t,  -/?o)  D 

- 2^r - -  Z, - 777 - Bln 


m= 1 


71—  1 

only  for  T  =  0 


4bko 


(42) 


The  other  elements  of  scattering  matrix  can  be  obtained  by  changing 
the  port  of  initial  excitation. 


3.  NUMERICAL  RESULTS 

The  proposed  Fourier  transform  technique  has  been  applied  to  the 
analysis  of  a  four-port  E-/H-plane  junction  and  a  three-port  E-/H- 
plane  junction.  The  results  have  been  compared  with  available  nu¬ 
merical  and  experimental  data.  The  scattering  parameters  were  cal¬ 
culated  from  (39)— (42).  The  modal  expansions  of  scattered  fields  in 
the  arm  waveguides  H  and  HI  were  truncated  at  m  =  n  =  M .  Then 
(27),  (30),  (33),  and  (34)  are  rendered  into  a  matrix  equation  with 
4 M(M  +  1)  unknowns.  We  first  consider  the  Magic-T  which  is  a  junc¬ 
tion  of  four  identical  waveguides  comprising  an  E-/H-plane  T-junction. 
The  values  of  geometrical  parameters  in  Figure  1  were  chosen  to  be 
2 a  =  2 w  =  15.799  mm  and  26  =  2t  =  d  =  c  =  7.899  mm.  Table  I 
shows  the  convergence  of  the  scattering  parameters  for  /  =  16.5  GHz 
and  the  CPU  time  of  computation  versus  the  truncation  number  M 
of  modal  expansions  in  waveguides  II  and  IE  .  The  numerical  compu¬ 
tation  was  run  on  Sun  Ultra  1  (UltraSPARC  143  MHz)  workstation. 
We  can  see  that  a  very  good  convergence  is  achieved  for  all  scattering 
parameters  as  the  number  of  modes  increases.  When  the  mode  number 
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M 

511 

Sn 

531 

5«: 

CPU(s) 

1 

0.097915!  23.79° 

0.5599531  -22.89° 

0.4191741  -109.64° 

0.7079251 17.66° 

0.02 

2 

0.159804161.49° 

0.5658491  -30.87° 

0.4367391-104.27° 

0.6808351 8.00° 

0.04 

3 

0.160002160.56° 

0.5663541  -30.96° 

0.4342701-104.17° 

0.6819471  7.83° 

0.12 

4 

0.163797162.27° 

0.563544!  -31.26° 

0.4352971-103.77° 

0.68271917.21° 

0.33 

5 

0.164247161.83° 

0.5636751  -31.28° 

0.4333201-103.86° 

0.68375917.18° 

0.79 

6 

0.165319162.44° 

0.5624501  -31.37° 

0.4336501-103.73° 

0.68430116.99° 

1.62 

7 

0.1655651 62.23° 

0.562507]  -31.38° 

0.4326341-103.79° 

0.6848371 6.98° 

3.50 

8 

0.166167162.49° 

0.5620271  -31.41° 

0.4327881-103.73° 

0.6849891 6.88° 

6.62 

9 

0.166295162.42° 

0.5620581  -31.42° 

0.4324091-103.75° 

0.6851711 6.88° 

12.70 

10 

0.166698162.56° 

0.5618451  -31.44° 

0.4324961-103.71° 

0.68519316.82° 

21.69 

11 

0.166787162.52° 

0.561864!  -31.44° 

0.4322581-103.72° 

0.6853061 6.82° 

37.73 

12 

0.167047162.61° 

0.5617091-31.45° 

0.4323121-103.70° 

0.6853351 6.79° 

62.49 

13 

0.167112162.58° 

0.5617221-31.46° 

0.4321211-103.71° 

0.6854301 6.79° 

97.46 

14 

0.167281162.64° 

0.5616011-31.46° 

0.4321581-103.69° 

0.68546416.76° 

153.41 

15 

0.167330162.62° 

0.5616101-31.47° 

0.4320101  -103.70° 

0.6855381  6.76° 

237.35 

Table  I.  Convergence  of  the  magnitude  and  phase  of  scattering  pa¬ 
rameters  at  /  =  16.5  GHz  versus  the  truncation  number  M  of  modal 
expansion  for  a  Magic  T-junction  with  2 a  =  2w  =  15.799  mm  and 
26  =  2t  =  c  =  d  =  7.899  mm. 


is  truncated  at  M  =  4,  the  relative  errors  of  the  solutions  are  about 
2%  for  1 5*11 1  and  less  than  1%  for  IS21I,  1*5311)  and  |S4i|.  It  is  worth 
noting  that  the  results  given  in  Table  I  satisfy  accurately  the  energy 
conservation  relation  ^t1=1  |*5u|2  =  1  for  any  truncation  of  modal  ex¬ 
pansion.  Figure  2  illustrates  |<Sii|,  |*52i|,  IS31I,  and  | *S*4i  |  as  functions 
of  frequency  for  different  numbers  of  modes  used  in  the  analysis.  The 
scattering  parameters  converge  very  fast  and  no  significant  difference 
is  observed  between  the  results  for  M  =  4  and  for  M  =  6.  Note  that 
the  total  number  of  modes  used  is  80  when  M  =  4.  The  same  struc¬ 
ture  has  been  analyzed  [3]  by  using  the  generalized  admittance  matrix 
approach.  For  the  sake  of  comparison,  the  numerical  results  of  |>Sii| 
given  in  [3]  are  reproduced  in  Figure  3.  The  curves  indicated  by  the 
labels  a,  6,  c,  d ,  e,  and  /  were  obtained  by  using  19,  24,  28,  38, 
45,  and  60  modes  in  each  of  four  ports,  respectively.  It  follows  that  the 
generalized  admittance  matrix  approach  requires  the  total  180  modes 
at  least  to  obtain  the  stable  and  convergent  results.  Comparing  Figure 
2(a)  with  Figure  3,  it  is  seen  that  the  total  mode  numbers  used  in 
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Figure  2d 

Figure  2.  Magnitudes  of  the  scattering  parameters  for  a  Magic-T 
junction  with  2a  =  2w  =  15.799  mm  and  26  =  2t  =  c  =  d  = 
7.899  mm  for  the  different  truncation  number  M  of  modal  expansion; 
(a)  | tS'n | ,  (b)  | *S,2i | ,  (c)  1 531 1,  (d)  |<S4i|.  The  total  number  of  modes 
used  in  the  analysis  is  4 M(M  +  1)  for  each  of  M. 
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Frequency  (GHz) 

Figure  3.  Scattering  parameter  |<Sii|  calculated  using  the  generalized 
admittance  matrix  approach  for  the  same  Magic-T  junction  as  shown  in 
Figure  2.  The  results  have  been  reproduced  from  Fig.  12(a)  of  reference 
[3].  The  number  of  modes  used  in  each  of  four  ports  are  19,  24,  28,  38, 
45,  and  60  for  the  curves  a,  b,  c,  d,  e,  and  /,  respectively. 

the  present  method  is  far  less  than  that  of  the  generalized  admittance 
matrix  approach  to  obtain  the  same  accuracy.  Another  important  dif¬ 
ference  is  in  the  relative  convergence  problem.  For  the  generalized 
admittance  matrix  approach,  there  exists  an  optimal  ratio  [3]  between 
the  number  of  TE  and  TM  modes  used  to  achieve  a  faster  and  stable 
convergence.  We  have  never  observed  such  a  situation  for  the  present 
method.  Using  the  Fourier  integral  representation  of  fields  in  the  main 
waveguide  I,  the  modal  expansion  coefficients  in  two  arm-waveguides 
H  and  IE  have  been  directly  related  through  the  algebraic  equations 
(27),  (30),  (33),  and  (34).  This  reduces  the  number  of  unknowns  to  be 
determined  and  removes  the  problem  of  relative  convergence. 

Figure  4  shows  the  numerical  results  of  the  magnitudes  of  the  scat¬ 
tering  parameters  for  the  Magic-T  with  2 a  =  2w  =  22.86  mm ,  2b  = 
2 1  =  10.16  mm,  and  c  =  d  =  11.43  mm.  The  results  are  compared 
with  the  measured  data  [2],  For  the  analysis  of  |<S22|  and  IS33I  in 
Figure  4(b),  we  have  changed  the  port  of  initial  excitation.  It  is  seen 
that  the  present  results  agree  very  well  with  the  measured  one. 

We  consider  next  an  asymmetric  E-/H-plane  three-port  junction. 
The  magnitudes  of  scattering  parameters  of  the  junction  with  2a  = 
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Figure  4.  Magnitude  of  the  scattering  parameters  for  a  Magic-T 
junction  with  2 a  =  2 w  =  22.86  mm,  2b  =  2t  =  10.16  mm,  and 
c  =  d  =  11.43  mm. 
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Figure  5.  Magnitude  of  the  scattering  parameters  for  an  asymmetric 
E-/H-plane  three-port  junction  with  2 a  =  2 w  =  15.799  mm,  2b  = 
2t  =  2d  =  c  =  7.899  mm.  The  left  side  end  of  waveguide  I  in  Figure  1 
is  short-circuited  at  z  =  0. 


2 w  =  15.799  mm  and  2b  =  2t  =  2d  =  c  =  7.899  mm  are  plotted  in 
Figure  5  for  the  excitation  of  each  of  three  ports,  and  are  compared 
with  those  [3]  obtained  by  the  generalized  admittance  matrix  approach 
using  28  modes  in  each  of  ports.  We  note  that  both  results  are  in  good 
agreement.  A  small  discrepancy  is  observed  in  | Sj  1 1  and  |S2i|,  since 
the  number  of  modes  used  in  [3]  was  not  sufficient  for  the  accurate 
computation. 

4.  CONCLUSION 

A  rigorous  and  efficient  technique  for  the  analysis  of  a  rectangular 
waveguide  E-/H-plane  junction  has  been  presented.  The  approach 
is  a  combination  of  the  Fourier  transform  technique  and  the  mode¬ 
matching  method.  By  introducing  an  idea  of  image  waveguide,  the 
problem  of  the  junction  was  transformed  into  a  problem  of  rectangular 
apertures  located  on  the  E- /H-planes  of  the  main  waveguide  of  infinite 
extent.  When  the  Fourier  transform  technique  is  applied  to  the  fields 
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of  the  main  waveguide,  a  set  of  linear  equations  relating  the  amplitudes 
of  the  normal  mode  expansion  in  each  of  arm  waveguides  are  obtained 
in  analytically  closed  form.  This  technique  reduces  significantly  the 
number  of  unknowns  to  be  determined  and  removes  the  problem  of 
relative  convergence.  The  scattering  parameters  can  be  calculated  by 
solving  a  matrix  equation  of  relatively  small  dimension.  The  numerical 
results  of  the  scattering  parameters  calculated  for  the  Magic-T  junc¬ 
tion  and  an  asymmetric  E-/H-plane  three-port  junction  were  in  a  very 
good  agreement  with  the  available  numerical  and  experimental  data. 


APPENDIX 

Using  the  residue-calculus,  the  integral  in  (28)  is  evaluated  as  follows: 


t(  1  +  5n0)5nq  .y,  Cpu  [1  ~  (-l)»+?e-^]  Qqn 

k0kpn  tan(2 kpnb)  y=0  k0b(l  +  5V o)  (C)L  ~  ?n)  ( Qv  ~  0 

—  1  _|_  g2 jCpvt  —  ^ —  (  —  ]_)^j  C~~ 

2bk0{l  +  8uo)(Cpu  ~  i2n)(Cpu  -  t\) 

(Al) 

with 


Qqn  — 


where  kpn  =  k$  —  cip  —  b^  and 
in  (Al)  is  calculated  from  the  residues  at  C  =  =t tn,  the  second  term  is 
from  the  residues  at  C  =  ±Cpu  which  satisfy  tan(2£p6)  =  0,  and  the 
third  term  is  from  the  residues  at  £  =  Cpu  ■  It  is  easy  to  show  that  the 
poles  of  the  integrand  in  (28)  located  at  (  =  ±tn  and  at  (  =  —(pu 
have  no  contribution  to  the  result  of  the  integration  for  a  finite  value 


A 

q  +  n  =  even 
A  .  . 
q  +  n  =  odd 


(A2) 


Cpu  =  \J kQ  —  ap  —  b2u.  The  first  term 
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of  r  .  Similarly,  the  integral  in  (29)  is  evaluated  as  follows: 
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(A3) 

where  ( pn  =  \J -a2-  and  gmn  =  \Jkl-w2m-  b„.  The  first 
and  second  terms  in  (A3)  are  calculated  from  the  residues  at  (  =  ±(pn 
and  the  third  term  is  from  the  residues  at  £  =  ±wm-  The  poles  of  the 
integrand  in  (29)  at  £  =  ±tn  has  no  contribution  to  the  result  of  the 
integration.  The  integrals  in  1%,  I4,  and  J\  to  J4  defined  by  (31), 
(32),  and  (35)  to  (38)  are  evaluated  in  closed  form  using  the  similar 
residue-calculus . 
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